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I .  Introduction 

We  consider  here  two  related  prohlems.  The  first  can  he  viewed  as 
an  extension  of  the  classic  two-sample  problem  which  asks  whether  two 
samples  have  the  same  underlying  distribution.  Thus,  given  two  indepen¬ 
dent  samples  . . . ,  . . . ,  where  the  X^  are  mutually  inde¬ 
pendent  with  common  c.d.f  F(x),  and  the  are  mutually  Independent  with 

common  c.d.f  G(x),  the  classic  two-sample  hypothesis  is  that  F(x)  equals 
G(x)  and  the  alternative  is  that  F(x)  and  G(x)  are  not  identically 
equal.  If  the  two-sample  hypothesis  is  true  and  (AX^  +  B)  is  substituted 
for  X^  (l  =  1, ...n)  (A  >  0,  B  real)  before  the  test  is  made,  any  consis¬ 
tent  test  for  the  classic  problem  would  reject  the  hj-pothesis  with  a  high 
probability.  The  numbers  A  and  B  represent  nuisance  scale  and  loca¬ 
tion  parameters  and  do  not  affect  the  form  of  the  distribution  of  (AX^  +  B) . 
Under  some  circumstances,  the  relevant  question  miglit  be  whether  F(x)  and 
G(x)  are  of  the  same  form  even  thou^  they  differ  through  the  presence  of 
the  above  mentioned  scale  and  location  paraimeters.  Formally,  the  hypothesis 
would  be  that  F(x)  =  G(Ax  +  B)  for  some  unspecified  pair  A,  B  (A  >  0,  B 
real).  The  alternative  is  that  the  above  equality  does  not  hold.  We  pro¬ 
pose  a  family  of  test  statistics,  all  members  of  which  lead  to  a  consis¬ 
tent  test  of  the  last  named  hypothesis  xinder  some  mild  restrictions  on  the 
form  of  F(x)  and  G(x) . 


Having  described  the  general  hypothesis  which  we  are  testing,  we 
shall  suggest  some  specific  instances  where  it  would  he  useful.  The 
range  of  application  will  he  seen  to  he  essentially  the  same  as  that  of 
the  classic  two-sample  test.  If  the  experimenter  wants  to  test  whether 
two  distrihutions  are  identical  and  he  knows  that  the  measurements  on  the 
two  populations  have  been  made  with  instruments  whose  zero  and  scale  cali¬ 
brations  are  different,  then  he  can  use  the  proposed  test  to  eliminate  the 
effect  of  these  extraneous  factors.  If  he  knows  the  numerical  values  of 
the  zero  and  scale  factors  and  can  adjust  the  data  to  a  common  zero  and 
scale,  then  the  usual  tests  would  be  more  efficient. 

The  experimenter  might  consider  two  populations  the  same  if  they 
both  have  distributions  of  the  siame  form  even  thougli  they  differed  in  a 
scale  and  location  parameter.  Here  again  the  proposed  test  is  applicable. 
For  example,  two  machines  mi^t  be  considered  interchangeable  if  some 
measurable  characteristic  of  their  output  has  the  same  distribution  for 
each  machine.  It  migb.t  not  be  feasible  however  at  the  time  of  making  the 
test  to  adjust  the  machines  so  that  the  scale  and  location  parameters  of 
the  output  distributions  are  the  same.  Alternatively,  through  some 
accident,  the  machines  might  have  been  adjusted  differently  at  the  time 
the  data  were  collected,  and  re-adJusted  subsequently  so  that  at  the  time 
the  statistician  gets  the  data  there  is  no  way  of  knowing  what  the  rela¬ 
tive  scale  and  location  settings  were  at  the  time  of  data  collection. 

Using  the  proposed  test,  such  data  would  not  have  to  be  discarded. 

A  further  use  of  the  test  would  be  to  detect  a  "linear  drift". 
Suppose  it  is  believed  that  observations  at  time  zero  have  a  distribution 
F(x)  (unspecified)  and  that  observations  at  some  fixed  later  time  can  be 


2 


considered  to  "be  of  the  form  AY  +  B  where  again  Y  has  the  distribution 
F(x).  Then  the  proposed  test  will  reject  the  hypothesis  if  either  the 
distribution  of  Y  is  not  F(x)  or  if  the  linear  model  AY  +  B  is  not 
correct . 

The  second  problem  is  that  of  estimating  the  scale  parameter  A, 
assuming  that  F(x)  =  G(Ax  +  B)  {A  >  0,  B  real).  Again  a  family  of  esti¬ 
mates  is  proposed,  each  member  of  which  provides  a  consistent  estimate  of 
A  under  the  same  mild  restrictions  on  F(x)  that  are  made  for  the  test 
of  hypothesis. 

This  estimation  problem  could  arise  in  several  ways.  After  accepting 
the  experimenter  might  want  an  idea  of  the  size  of  the  scale  parameter 
A.  Alternatively,  ass\imlng  to  be  true  on  the  basis  of  theoretical 

considerations,  "A"  might  represent  some  physical  characteristic  which 
the  experimenter  wants  to  measure.  Tf  he  is  unwilling  to  make  specific 
assumptions  about  the  type  of  distribution  involved,  the  proposed  esti¬ 
mates  with  their  fairly  v;eak  assumptions  can  be  used.  If,  alternatively, 
the  experimenter  does  know  more  about  the  distribu lions,  br  ■  "inds  compu¬ 
ting  "best"  estimates  is  too  difficult,  the  proposed  estimates  might 
again  be  useful. 

A  common  problem  where  such  an  estimate  is  needed  is  in  the  case 

where  a  time  factor  introduces  a  linear  drift  into  the  observations.  If 

all  the  X's  are  assumed  to  have  been  measured  at  time  t  and  the  Y's 

o 

at  another  time  t^,  then  the  assumed  model  would  apply. 

In  Section  2,  we  describe  the  proposed  statistics  in  detail  and 
discuss  their  properties.  Section  5  consists  of  an  example  of  the 
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computations  Involved.  The  following  Section  is  concerned  chiefly  with 
establishing  the  consistency  of  the  proposed  estimates  and  tests.  The 
important  OLuestion  of  the  distributions  or  even  the  asymptotic  distri¬ 
butions  of  the  proposed  statistics  is  discussed  only  briefly  in  Section 
5.  Section  6  raises  some  additional  questions. 


k 


2. 


The  Test  and  The  Estimator 


Let  X^, y^  be  two  independent  sets  of  mutually 
independent  random  variables.  Let  F(x)  be  the  common  c.d.f  of  the 
X^  (i  =  l,...,n)  and  G(x)  be  the  common  c.d.f  of  the  y^  (i  =  l,...,n). 
We  shall  list  in  Section  4  the  necessary  regularity  conditions  on  F(x) 
and  G(x).  Let  XJ  <  Xi  <  . ♦ •  <  X'  be  the  ordered  values  of  Xt,...,X 
and  y|  <  y^  <  . . .  <  y^  be  the  ordered  values  of  y^, ...,  y^.  Define 
as  (X|^^  -  X|)  and  as  -  y^  (l  =  1,...,  n-l).  Let  a, 

Pj  V  be  arbitrarily  chosen  numbers  satisfying  the  following  res¬ 
trictions: 

(2.1)  0<a,  p<l}  0<B<i>  0<u<v<l  . 


We  then  define  the  following  statistics: 

*  *'1 

(sin  ttp/n«p)  2- 

l=[n^/2  t  6^ 


(2.2)  (p,6)  =  < 


w  1/2  +  6i 

i=[n-n  '  ] 

(l/n  log  n)  1. 

l.tnVs  .  6, 


^  [nv-1 ] 

(sin  itp/njtp(v-u))  ^ 

i= [nu+1  ] 


(2.5)  Sjj(p,u,v)  =  / 


[nv-1  ] 

(l/n  (v-u)log  n)  ^ 
i=[nu+l] 


(U^/V^)P  0<p<l 


(U^/V^)  P  =  1 


(U^/Vj)^  0<p<l 


(U^/V^)  P  =  1 
where  W  is  the  greatest  integer  less  than  or  equal  to  Ti , 
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similarly,  define  S'(p,6)  and  S'(p,u,v)  'byreplacing  (ll/V.)  by 

n  n  11 

(■y^/U^)  In  the  corresponding  definitions  for  S^(p,6)  and  S^(p,u,v). 

We  desire  to  test  the  hypothesis 

H  :  F(x)  =  G(Ax  +  B) 

(2.4) 

:  F(-x)  ^  G(Ax  +  B). 

The  first  proposed  family  of  tests  is  to  reject  If  the  product 

Sjj(p>8)S^(p,8)  is  "too  large".  Cie  second  proposed  family  of  tests  is 
to  reject  if  the  product  S^(p,u,v)S^(p,u, v)  is  "too  large".  When 

F(x)  and  G(x)  satisfy  conditions  (4.2)  to  (4.4),  we  show  in  Section  4 
that  S^(p,8)  S^(p,6)  converges  stochastically  to  unity  when  is 

true  and  to  a  quantity  greater  than  unity  otherwise.  This  implies  that 
the  proposed  tests  are  consistent  against  alternatives  F(x),  G(x)  satis¬ 
fying  these  conditions.  We  show  also  that  the  tests  are  consistent 
against  distributions  which  satisfy  only  condition  (4.2).  Derivations 
from  condition  (4.2)  can  lead  to  situations  in  which  is  not  true 
and  the  probability  of  rejecting  will  not  tend  to  unity  as  n 
Increases.  This  can  happen  when  one  of  the  distributions  has  a  flat 
section.  An  example  is  given  by  (4.50).  If  is  true  and  conditions 
(4.2)  to  (4,4)  are  violated,  the  test  will  tend  to  accept  in  some 
cases  and  reject  in  others. 

Under  the  weaker  assumptions  (4.1)  the  proposed  tests  based  on 
Sjj(P;U,v)  S^(p,u,v)  can  be  shown  to  be  consistent  by  the  same  type  of 
argument  as  above.  Again,  violation  of  this  condition  will  lead  to  the 
situation  described  above.  Although  the  S^(p,u,v)  family  of  tests 
requires  weaker  assumptions  regarding  the  form  of  F'(x)  than  does  the 
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S^(p>B)  family  to  guarantee  the  same  operating  properties,  the 
S^(p,u,v)  teats  require  that  a  fixed,  percentage  of  the  observations  must 
be  Ignored,  whereas  the  percentage  ignored  by  the  Sj^(p,B)  family  goes 
to  zero  with  n.  Both  families  of  tests  are  invariant  under  separate 
linear  transformation  on  the  X's  and  Y's.  This  is  a  reasonable  require¬ 
ment  for  any  test  of 

The  choice  of  p  would  presumably  be  dictated  by  power  considera¬ 
tions.  Power  depends  on  distributions,  and  we  know  almost  nothing  about 
these. 

Assuming  that  F(x)  =  G(Ax  +  B),  we  propose  to  estimate  the  scale 
factor  A  by  S^(p,S)  if  assumptions  (4.2)  to  (4,4)  are  satisfied,  or 
4y  S^(p,u,v)  if  condition  (4.1)  is  satisfied.  The  consistency  of 
these  estimates  follows  from  Theorem  2,  Lemma  2,  and  Theorem  3.  Again 
choice  of  p  depends  on  the  unknown  distributions  of  the  statistics. 
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3.  Computational  Example 

To  Illustrate  the  type  of  computations  involved,  we  shall  compute  the 
scale  parameter  estimate  S^(p>  •!>  *9)/  with  n  =  10  for  convenience. 

The  X's  represent  ten  observations  which  were  taken  from  a  uniform  dis¬ 
tribution  on  the  unit  interval  by  means  of  a  table  of  random  numbers. 

The  Y's  represent  ten  observations  from  a  uniform  distribution  on  the 
Interval  from  zero  to  five.  Thus  the  factor  A  to  be  estimated  is  five. 
The  table  below  shows  the  steps  needed  to  go  from  the  and  to  the 
Xj^  and  (ordered  values)  to  the  and  (differences)  to  obtain 
finally  the  estimator  S^(p,u,v).  Two  values  of  p  are  used,  p  =  l/2 
and  p  =  1  to  illustrate  the  two  types  of  estimators  described. 


i 

"1 

^i 

^1 

Y' 

i 

'^i 

V^i 

1 

.725QC 

.05275 

.13538 

4.44515 

0.06115 

0.42773 

3.160 

1.778 

2 

.30061 

.18813 

.04203 

0.72285 

0.48890 

0.2339; 

5.566 

2.559 

5 

. 31694 

.23016 

.05608 

2.78165 

0.72285 

0.62230 

11.097 

3.331 

4 

.18813 

.28624 

.01437 

4.51455 

1.34515 

1.43650 

99.965 

9.998 

5 

.05275 

.30061 

.01633 

0.06115 

2.78165 

0.57620 

55.285 

5.940 

6 

.79607 

.31694 

.03206 

1 

4.77130 

3.35785 

1.08730 

33.914 

5.824 

7 

.34900 

.34900 

.039^ 

0.48890  ! 

4.44515 

0.06940 

1.761 

1.327 

8 

.23016 

.38840 

•33746 

4.82833 

4.51455 

0.25675 

0.761 

0.872 

9 

.38840 

.72586 

.07021 

1.3451:' 

4.77130 

0.05705 

0.813 

0.902 

10 

.28624 

.79607 

~ 

5.35785 

4.82835 

— 

— 

— 

From  (2.3),  we  see  that  to  compute  •9)>  we  must  sum 


(Vg/Ug)  through  (Vg/Ug),  and  divide  the  total  by  (.8)(l0)  log  10,  which 
is  I8.t21.  The  total  is  188.350  and  thus  our  estimate  S|q(1,  .1,  .9) 
is  10.225.  To  obtain  S^q(1/2,  .1,  .9),  we  sum  (Vg/Ug)^/^  throng 
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and  divide  "by  (.£)(lO)(rt/2)  (since  sin  rt/2  =  l)  which  is 
12.566.  The  total  here  is  29.65I  and  thus  our  estimate  Sj^Q(3./2,  .1,  .9) 
is  2.360.  It  is  seen  that  one  estimate  is  more  than  double  the  true 
value  while  the  other  is  less  than  half.  Thus,  neither  is  particularly 
good.  Since  the  seunple  size  is  quite  small,  we  would  not  expect  the 
asymptotic  properties  of  the  estimator  to  have  any  bearing  and  this  ex¬ 
ample  can  be  considered  as  an  illustration  of  the  danger  of  using  a  small 
sample  when  the  small  sample  properties  of  the  statistic  are  not  known. 
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V.  Convergence  of  the  StatleticB 


Let  X,,...,X  }  Y.,...,Y  and  X*  <  ...  <  X'  ;  and  Y'  <  ...  <  Y' 

1-  -n'  1-  -n 

be  as  described  in  Section  2.  We  shall  list  below  regularity  conditions 
on  the  distributions  F(x)  and  G(x).  These  conditions  will  sometimes 
overlap,  and  not  all  of  them  need  be  satisfied  for  the  stochastic  con¬ 


vergence  of  a  particular  family  of  statistics.  We  state  the  conditions 
in  terms  of  F(x)  for  convenience. 


(4.1)  For  a  given  pair  of  numbers,  u  ,  v  (0<  u  <v  <  l)  F’^(u  )  and 

0  0—0  o—  o 

F"^(v^)  are  uniquely  determined,  and  F(x)  has  on  [F“^(u^, 
F'^(v^)]  a  derivative  f(x)  \diich  has  only  a  finite  number  of 
discontinuities,  is  bounded  above,  and  is  bounded  away  from  zero. 


(4.2)  F(x)  has  a  derivative  f(x)  which  is  bounded  for  all  x,  and 
has  only  a  finite  number  of  discontinuities.  Further,  there  ex¬ 
ists  a  positive  value  D'  such  that  for  any  value  D  in  the 
open  interval  (0,  D')^  the  set  of  points  (x  ;  f(x)  >  D)  is 
an  interval. 

(4.5)  If  f[F”^(r)]  is  not  bounded  away  from  zero  as  r  approaches 
unity  then  either  there  is  a  value  t(0<  t  <l)  such  that  the 
quantity  d  f[F“^(r)]/dr  is  negative  and  non- increasing  for 
r  >  t,  or  the  quantity  d  f[F“^(r)]/dr  approaches  a  finite  ne¬ 
gative  limit  as  r  approaches  unity. 

(4.4)  If  f[F"^(r)]  is  bounded  away  from  zero  as  r  approaches  zero 

then  either  there  is  a  value  t'(0<t'<l)  such  that  the  quantity 
d  f  [F''^(r)  ]/dr  is  positive  and  non-decreasing  for  r  <  t',  or 
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the  quantity  d  f[F"^(r)/dr  approaches  a  finite  positive  limit 
as  r  approaches  zero. 

Note  that  implies  (^.l).  Also,  the  second  part  of  (4.2)  is 

satisfied  by  every  unimodal  f(x).  Conditions  (4.5)  and  (4.4)  can  be 
verified  by  using  the  fact  that  d  f[F"^(r)]/dr  =  d  l(ig  f [F"^(r) ]/dF"^(r) . 
Thus  the  condition  (4.3)  says  that  if  f(x)  goes  to  zero  as  x  increases, 
then  there  exists  a  finite  y  such  that  d  log  f(x)/dx  is  non- increasing 
for  X  >  y.  The  conditions  (4.3),  (4.4)  are  satisfied  by  all  infinite 
Polya  frequency  functions.  The  exponential,  normal  and  Welbull  distri¬ 
butions  are  well-known  examples. 

Now  we  shall  consider  the  question  of  stochastic  convergence.  Let 
h(t)  be  a  bounded  non-negative  function  of  t  defined  for  0<  t  <  1. 

For  p  in  the  interval  define  H^j(p)  follows: 

H_(p)  =  (1/n)  h(i/n)(U /V  )P  ,  0<  p<  1 

1=1  ^ 

(4.5) 

n-1 

H  (p)  =  (l/n  log  n)  Y.  4(i/n)(u7V. )  p  =  1. 

i=l 


Then  we  have 

Theorem  1.  If  f(x)  =  G(x)  =  x,  (0<  x<  l),  then  as  n  increases, 
Hj^(p)  converges  stochastically  to 

n-1 

(np/n  sin  np)  ^  h(i/n)  0<p<  1 

1=1 

(4.6)  n-1 

(l/n)  y  h(i/n)  p  =  1  . 

i=l 


I 
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Proof; 


Let  ^  and  Z^, ^  be  two  independent  sets  of 

*x 

mutually  independent  random  variables  with  cdf  e  .  Let  R  = 

n-1  i=l 

T  =  y  2  .  Then  it  is  well  known  that  (W,/R  ,  ...,W  -./R  )  and 

ti  ^  ^  *1  n*"  iL  n 

(U^, have  the  same  Joint  density.  Similarly  for 
(Z^/T^,  . .  .,Z^  ^^l'***’^n  1^*  Independent  of 

(Wi/R^,  . .  .,Wn_i/Rjj)  and  is  independent  of  • 

Thus  we  need  only  consider  the  convergence  of 


(^.7) 


n-1 


(1/n)  h(i/n)  (T^  W.^/Rj^  Z^)^ 


0  <  p  <  1. 


Using  the  strong  law  of  large  numbers ,  and  making  somo  straight¬ 
forward  computations,  we  find  that  converges  to  unity  w.p.l 

for  all  p,  0<  p  <  1.  Also,  for  0<  p  <  -^  ,  we  find  that 

(4.8)  (1/n)  h(i/n)  (W./Z,)5 

i=l  ^ 

converges  to  (4.6)  w.p.l.  For  ■^  <  p  <  1,  the  results  in  Gnedenko  and 
Kolmogorof  [2]  on  the  relative  stability  of  a  sum  of  positive  random 
variables  with  finite  expectations  show  that  (4.8)  converges  stochasti¬ 
cally  to  (4.6). 

For  the  case  p  =  1,  we  must  consider 

n-1 

(4.9)  (l/n  log  n)  Y.  ^(i/n)  (WVZ. )  . 

i=l  ^  ^ 

Again  the  convergence  of  (4.9)  stochastically  to  (4.6)  follows  from 
results  in  [2]  on  asymptotic  stability  of  sums  of  positive  random 
variables.  The  boundedness  of  h(t)  is  important  in  all  of  the  above 
convergence  computations.  This  proves  Theorem  1. 
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Corollary  1.  If  P(x)  =  G(x)  =  x  (0<x<l),  and  I  h(t)dt  exists 

”  ”  •'o 

(in  the  Eiemann  sense),  as  n  Increases,  converges  stochastically 


to 


(U:.10) 


r 


(np/sin  itp)  J  h(t)dt 


I 


h(t)dt 


0  <  p  <  1 


P  =  1 


Proof: 


If  I  h(t)dt  exists,  then  as  n  Increases,  approaches  (4.10). 

•^o 

Hiis  proves  the  corollary. 

Nov  suppose  condition  (4.1)  is  satisfied  and  let  8^(p,u,v)  be  de¬ 
fined  by  (2.3)  with  0  <  <  u  <  v  <  <  1,  where  u^,  v^  are  given 

in  (4.1)  Then  we  have 

Hieorem  2.  As  n  increases,  S  (p,u,v)  converges  stochastically  to 


(4.11) 


(l/(v.u))J^  (g[G"^(t)]/f[F"^(t)])P  dt 


Proof: 

If  and  X^  are  both  in  [F"^(u),  F"^(v)],  and  and  Y^ 

are  both  in  [g“^(u),  G“^(v)],  we  can  write 

F(X',i)  -  F(x-)  =  f(0j)(X-,,  -  Xj)  ,  (X«  <  Oj  <  X',,) 

(4.12) 

-  jj)  ,  (jj  <  • 
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Define  \  as  |P(x;)  -  (j/n)l  and  K  as 
<J  J  J 


|G(Y^)  -  (d/n) I  .  By  the 

i/a-6. 


Gllvenko-Cantelli  theorem,  we  know  that  for  any  positive  B,  max  n' 

i 

converges  stochastically  to  zero  as  n  increases.  Define  6.  as 

J 

|X^  -  F'^(j/n)| .  If  Xj  and  F“^(j/n)  are  both  in  the  interval 
[F  ^(u),  F“^(v)],  then  since  f(x)  >  A  >  0  on  this  interval  (by  (k.l)), 


we  have  6j<(Xj/a).  Then  if  X',  F‘^(j/n)  and  F"^( (j+l/n)  are 

all  in  the  interval  [F"^(u),  F"\v)],  we  have  that  I©,-  F"^(j/n)l  <  8,+ 

d  *"  d 


+  l/nA,  and  we  can  wite 


(4.13)  F(X^^^)  -  F(X')  =  f[F-^(j/n)]  (Xj^^  -  Xp  +  -  X') 


where  =  f(0j)  -  f[F“^(j/n)).  But  because  of  the  uniform  continuity  of 

f(x)  in  fF  ^(u),  F  ^(v)],  the  above  bound  for  |©j  -  F“^(j/n)|  and  the 

Glivenko-Cantelli  theorem,  it  is  easily  seen  that  max  |y. |  converges 

nu<  i  <  nv 

stochastically  to  zero  as  n  increases.  Similarly  we  can  show  that 


(4.14)  ®(^j+l)  -  =  g[G"^(d/n)](Y^^^-  Yp  + 


where  max  ly'l  converges  stochastically  to  zero  as  n  increases, 
nu  <  i  <  nv 

Denote  F(X^^^)  -  F(Xp  by  Wj(j  =  l,...,n-l),  and  -  G(Yp 

(J  =  l,...,n-l)  by  Zj.  Notice  that  (W^, . . .,W^_^)  and  (Z^,...,Z^  ^) 
are  distributed  as  two  independent  sets  of  sample  successive  differences 
from  the  uniform  distribution  on  the  unit  interval.  Using  (4.13)  and 
(4.14)  we  have 

(^•15)  I  (g[G“^(i/n)]/f[F-l(i/n)])P  (W /Z  )5  =  E  (U,/V,  )P  + 

nu<i<nv  ^  ^  nu<i<nv  ^ 

+  E  (U,/V  )^  {(f(0  )  g[G■^(i^)]/g(0')f[F■^(i/n)|)P  -1)  . 
nu<  i  <  nv  ^  ^  ^  J 
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Using  the  result  of  Corollary  1,  we  have  that  the  left  side  of  (4.15) 
when  properly  normalized  converges  stochastically  to  (4.11).  Let  the 
appropriately  normalized  second  term  on  the  right  side  of  (4.15)  he  re¬ 
presented  hy  S^(p,u,v).  As  n  increases,  the  probability  approaches 
\inity  that  ls^(p,u,v)|  will  be  no  greater  than 

S  (p,u,v)  {(  max  +  max  |7*^|b^)/A^C^) 

nu<i  <  nv  nu<  i  <  nv 

where  0  <  A  <  f(x)  <  B  <  »  for  F“^(u)  <  x  <  F"^(v),  and  0  <  C  <  g(x)<  D  <  «, 

for  G"^(u)  <  X  <  g”^(v).  Ihus  we  have  that  |§jj(p,u,v) |/s^(p,u,v)  con¬ 
verges  stochastically  to  zero  as  n  increases.  But  S^(p,u,v)  +  S^(p,u,v) 
converges  stochastically  to  (4.11)  as  n  increases.  This  proves  Theorem 


2. 

It  is  clear  that  the  same  proof  will  show  that  under  condition  (4,l), 
S'(p,u,v)  (defined  in  Section  2)  converges  stochastically  to 


4.15)  (l/(v-u))  (f[F“^(t)]/g[G"^(t)])^  dt. 

-^u 


Thus,  under  condition  (4.l),  we  obtain  that  as  n  increases  the 
statistic  S^(p,u,v)  S^(p,u,v)  converges  stochastically  to 

V  V 

(4.16)  (iMv.u)f(f  (f[F"^(t)]/g[G"^(t)])Pdt)(f  (g[G”^t)]/f[F-^(t)])J’dt). 

-'u  “^u 

Nov,  assume  that  conditions  (4.2)  to  (4.4)  are  satisfied.  Let 
Sj^(p>B)  be  given  by  (2.2).  Then  we  have 

Theorem  ?.  If  F(x)  =  G(Ax  +  B),  then  as  n  increases,  S^(p,6)  con¬ 
verges  stochastically  to  (l/A)^  . 
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Proof! 


We  assume  that  either  or  both  of  the  quantities  11m,  f(x), 

xAp"^(0) 

lim  f(x)  is  zero.  Otherwise,  nieorem  2  applies.  For  arbitrarily 
xtF"^(l) 

chosen  u,  v,  0  <  u  <  v  <  1,  we  have 


(^■19)  S^(p,8)  =  (v-u)  S^(p,u,v)  +  b(n,p)  /  *  Z 


where 


(4-.ao)  b(n,p) 


sin  itp/nnp 
l/n  log  n 


nl/^^®<i<  nu  nv  <  i<  n-n^^'*‘®i 


0  <  p  <  1 

P  “  1 


By  Theorem  2,  (v-u)  Sj^(p,u,v)  converges  stochastically  to 

J  (g[G“^(t)]/f[P"^(t)])^  dt  -  (v-u)(l/A)^  . 


We  shall  now  investigate  (U^/V^)®  ,  the  treatment  of 

nv<  i<  n-n^/^"^® 

Z  (u^/v.)P 

jjV2+B  ^  being  entirely  similar. 


Since  F(x)  =  G(Ax  +  B),  we  see  that  has  the  same 

distribution  as  (AYj^  +  B, +  B)  where  ?  has  the  distribution 
F(x).  Letting  the  ordered  be  <  ...  <  ?^  ,  and  the  differences  be 
^i  we  have  that  ^  (Uj/V^)^  has  the  same 

distribution  as  J]  (U^/AV|)®  .  Thus  we  want  to  study  the  convergence  of 

y (u  /v« )P 

^  ^  ^  1/2+8  *  condition  (4.3)  and  the  relation  (4.12),  we  ob- 

nv<  i<  n-n  ' 

tain  that  y  (U^/v^)^  can  be  bounded  above  by 

nv<  i<  n-n^/®"^* 
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(^•22)  (f(Y;)/f(X'^^))P{(F(X^^^)(-P(Xp)/(P(?'^^)-p(?p))P 

If  we  can  find  a  bound  (say  B)  for  f (?p/f (X|^^),  then  we  would  be 
through,  for  (^.22)  woxild  then  be  bounded  by 

(4.25)  ®*’nv<^<n  * 


By  Corollary  1,  (4.23)  when  normalized  to  agree  with  S^(p,B)  con¬ 
verges  stochastically  to  (l-v)B^.  Thus,  by  choosing  v*  (say)  sufficiently 
large,  we  will  have  that 


approaches  unity  as  n  increases  for  each  v  such  that  v*  <  v  <  1. 

To  bound  f  (Xp/f(Xj^_j_j^),  we  note  that  If  >  then  unity 

serves  as  a  bound.  Thus  we  must  find  a  bound  (B')  for  this  ratio  when 
^1  ^  ^i+1  ‘  then  take  B  =  1  +  B' .  Letting  t^  P(?^)  and 

tg  =  F(X^^^)  (momentarily  suppressing  i),  we  have  t^  ^  tg  ,  and 

(4.24)  (f(5p/f(X|^^))  =  (f[P"^(t^)]/f[P'^(t2))). 

By  the  mean  value  theorem,  and  the  fact  that  f[P"^(l)]  »  0,  we  have 


f  [P"^(t^)]  =  (l-t^)(-df[P"^(t)]/dt)^^g  <  9l  <  1 


(4.25) 


f[P''^(tp)]  =  (1-t  )(-df[P"^(t)]/dt)._ 


t=02  ^2  5  ®2  ^  ^  • 


From  condition  (4.3),  we  see  that  there  Is  a  value  v^^  such  that  if 

*1*  ^2  exceed  v^,  then  either  (df (P“^(t) ]/dt)^^g  (j=l,2)  is 

J 


♦ 
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sufficiently  near  a  limit  L,  in  which  case  the  ratio  (4.8lf)  is  hounded 
by  2(l-t^)/(l-tg),  or  df [P"^(t)]/dt  is  a  decreasing  function  of  t  for 
t  >  t^.  In  the  latter  case,  by  Lenuna  1  below,  we  have  (l-t^)/(l-tg)  as 
a  bound  for  the  ratio  (4.2l^).  Clearly,  if  v  >  v^,  then  as  n  Increases, 

P  >  v^,  ,  nv  <  1  <  n} 

approaches  unity.  Thus,  we  need  only  to  find  a  bound  for 
[1-F(y^)]/[1-Pi[j£j^j^)  ],  (nv  <  1  <  n.-n^^'*’®).  ReVrltlng  this  last  ratio  as 

(i.r(5;)) 

n  n  1 

(a:^)  „VS-6  .. 

using  the  Glivenko-Cantelli  Theorem  (see  Theorem  2),  and  noting  that 
[(n-i)/n  '  ]  is  greater  than  unity,  we  have  that 

P  {([l-F(?^)]/[l.F(X^^j^)]  <  2  ,  nv  <  i  <  n-n^/^"^®} 

approaches  unity  as  n  increases. 

Putting  together  all  of  the  above  pieces,  and  noting  that  we  can 
choose  u,v  so  that  (v-u)  is  arbitrarily  close  to  unity,  we  see  that 

P  (|s^(p,6)  -  (1/A)^  I  <  €) 

approaches  unity  as  n  increases.  This  proves  Theorem  3. 

Lemma  1.  Suppose  f(x)  is  continuous  and  differentiable  on  the  Interval 
(a,b),  with  f'(x)  <  0  and  decreasing.  Let  <  Xg  be  in  (a,b).  By 
the  mean  value  theorem. 


f 
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(4.26)  f(x^)  =  f(b)  +  -  t)  f'(0^)  <  0^  <b  ,  i  =  1,2  . 

Then  ®2  -  ®1  * 

Proof; 

First  note  that 

(4.27)  f(x^)  =  fCxg)  -  (xg  -  x^)  (f'(05))  *1  <  <  Xg 

and  that  ®1  ®3  ' 

f(x^)  =  f(h)  +  (b.X3_)(-f'Cei))  <  f(fc)  +  (b-x^)(-f(05))  = 

(4.28) 

=  f(b)  +  (b-X2)(-f' (0^))  +  (Xg  -  x^)(-f (0^))  . 

Using  (4.27)  in  (4.28),  and  noting  (4.26)  we  have  0^  >  0g  which 
is  impossible.  Writing  f(x^)  =  f(kg)  +  (f(x^)  -f(Xg)),  then  using 

(4.26)  to  represent  f(xg)  and  (4.27)  to  represent  (f(x^)  -  f(x2))  and 
subtracting  the  result  from  the  representation  of  f(x^)  given  by 

(4.26) ,  we  obtain 

(4.29)  0  =  (x^  -  b)  f'(0j^)  -  (xg  -  b)  f'(a2)  +  (xg  -  x^)f'(0j) 

=  (xg  -  b)(f'(0^)  -  f'(e2))  +  (*2  - 

Since  9^  <  >  f'(9^)  is  positive,  or  9^  <  9^.  This  proves 

LemmEi  1. 

In  the  s€une  way  that  we  proved  'Hieorem  2,  we  can  show  that  under 
the  same  conditions,  S^(p,6)  converges  stochastically  to  A^.  Putting 
together  these  results,  we  have  that  Sj^(p,8)  S^(p>8)  converges  stochas¬ 
tically  to  unity  as  n  increases,  when  F(x)  =  G(Ax  +  B),  and  conditions 
(4.2)  to  (4.4)  are  valid. 
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Ve  shall  now  state  a  Lemma  due  to  Weiss  [3],  which  will  be  useful 


below. 

Lemma  2.  If  F(x)  and  G(x)  are  two  distribution  functions  and  u, 

V  (0  <  u  <  V  <  l)  are  two  given  numbers,  suppose  P"^(u),  P"^(v),  G'^tu), 
G"^(v)  are  all  uniquely  determined.  Also  suppose  that  F(x)  has  a 
derivative  f(x)  between  F“^(u)  and  F“^(v),  and  0(x)  has  a  deriva¬ 
tive  g(x)  between  G”^(u)  and  G"^(v).  Then  a  sufficient  condition 
that  f[F"''(r)]  =  kg[G“^(r)]  for  almost  all  r  in  [u,v]  (where  k  is 
a  positive  constant)  Is  that  there  are  two  constants  C,  D  (C>  0),  such 
that  r(Cx  +  D)  =  G(x)  for  all  x  in  the  interval  [G“^(u),  G'^Cv)]. 

If  in  addition,  f(x)  >  0  between  F"^(u)  and  P”^(v),  the  condition 
is  necessary  as  well  as  sufficient. 

We  omit  the  proof  since  it  is  contained  in  [3]»  Note  that  Weiss 
omitted  the  condition  f(x)  >  0  between  F"^(u)  and  P’^(v)  and  without 
this  condition  the  statement  is  incorrect.  A  simple  counter-example  is 
given  by  the  following  pair  of  distributions: 
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Now  we  shall  return  to  the  prohlem  of  examining  S^(p,&)  S^(p,8) 
when  conditions  (4,2)  to  (4.4)  hold  hut  F(x)  *  G(Ax  +  B).  We  have 
Theorem  4.  If  F(x)  *  G(Ax  +  B)  (A  >  O),  and  conditions  (4.2)  to  (4.4) 
obtain,  then  as  n  increases,  for  every  e  >  0  , 

(4.31)  P  (S^(p,S)  S^(p,6)  +  €  >  [  /  (frP"^(t)]/g[G"^(t)])^  dt]- 

r\ 


•if  (g[G”^(t)]/f[F"^(t)])^  dt]] 


approaches  unity. 

Proof! 

Let  S^(p,8)  and  S^(p,8)  he  represented  as  in  (4.19).  Then 
note  that  S^(p,5)  S^(p,8)  can  he  written  as 

(4.52)  (v-u)^  3^{p,u,v)  S^(p,u,v)  +  B(n,p,u,v) 

where  u,v  (0  <  u  <  v<l)  are  arbitrary  and  B(n,p,u,v)  is  strictly 
positive.  By  (4.17)  we  have  that  (v-u)^  S^(p,u,v)  S^(p,u,v)  converges 
stochastically  to 

,  V 


(4.35)  f  h(t)  dt  f  (l/h(t))  dt 

-^u  "^U 


where 


(4.34)  h(t)  =  (f[F"^(t)]/g[G"^(t)])^  . 


But  by  the  Lehesgue  Monotone  Convergence  theorem,  we  know  that  by 
making  (v-u)  close  to  unity,  we  can,  for  any  positive  e,  make  (4.33) 
greater  than 
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approaches  unity  as  n  increases.  This  establishes  consistency  when 
conditions  (4.2)  to  (4,4)  are  satisfied.  It  is  easily  seen  that  condi¬ 
tions  (4.5)  and  (4.4)  were  not  essential  to  the  argument.  Thus  even  if 
only  condition  (4.2)  holds,  the  test  will  be  consistent.  If  condition 
(4.2)  is  violated,  we  must  allow  distributions  such  as  those  given  in 
(4.30)  and  for  such  distributions  it  can  not  be  shown  that  the  probability 
(4.57)  approaches  unity  as  n  increases. 
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Remarks  About  Large  Sample  DlatrlTautlone . 


We  shall  give  a  heuristic  argument  to  indicate  vhy  for  some  values 
of  p,  S^(p,u,v)  and  S^(p,u,v)  might  have  a  bivariate  normal  limiting 
distribution,  whereas  for  other  values  of  p  this  distribution  is  not 
possible.  Similar  statements  can  be  made  about  S^(p,8)  and  S^(p,S). 
From  the  proof  of  Theorem  2,  we  see  that  S^(p,u,v)  and  S^(p,u,v)  have 
approximately  the  same  joint  distribution  as 

nv  .  . 

l'(njp)  {f[F"^(i/n)]/g[G"^(i/n)]}^  (W^^/Z^)^  and 

l=nv 

(5.1) 

1  i  n  r. 

b(n,p)  t  {g[G-^(i/n)]/f[F  (i/n)]}^  (Z^/W^)^ 

i=nu 

where  (W^, . . and  (Z^, . . . have  the  same  joint  distribu¬ 
tions  as  the  corresponding  quantities  in  Theorem  2,  and  b(n,p)  is 
given  by  (It’.SO) . 

From  the  proof  of  Theorem  1,  we  see  that  the  qviantities  in  (5*1) 
have  approximately  the  same  joint  distribution  as 

^(n,p)  ^  (h(i/n))^  and 

isnu 

(5-2) 

nv  „  _ 

b(n,p)  i  (l/h(l/n))P  (yX^)^ 

l=nu 

where 

(5.5)  h(i/n)  =  f[F-l(i/n)]/g[G"^(i/n)] 


I 
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and  X^)  are  Independent  random  variables  with  c.d.f  e“* 

and  (Y^, Y^)'  are  Independent  random  variables  with  o.d.f  e"^. 

Each  of  the  quantities  in  (5«2)  is  a  linear  combination  of  indepen¬ 
dent  random  variables,  and  when  p  <  l/2,  the  first  two  moments  of  these 
variables  are  finite.  Thus,  the  bivariate  central  limit  theorem  applies 
and  the  expressions  in  (5.2)  have  a  limiting  bivariate  normal  distribu¬ 
tion.  Thus,  for  p  <  1/2,  it  is  reasonable  to  suppose  that  S^(p,u,v) 
and  S^(p,u,v)  have  a  limiting  bivariate  normal  distribution.  When 
p  =  1/2,  all  moments  of  (X^/X^^)  up  to  the  second  exist,  and  each  term 
in  (5.2)  can  be  shown  to  have  a  limiting  normal  distribution.  Similarly, 
a  bivariate  limiting  normal  distribution  for  the  terms  in  (5.2)  can  be 
obtained  from  a  generalization  of  the  one-dimensional  result. 

When  p  >  1/2,  not  all  moments  of  (X^/X^)^  of  order  less  then 
two  exist.  If  one  of  the  expressions  in  (5*2)  had  a  limiting  normal 
distribution,  then  the  expression 

nv 

(5.4)  b(n,p)  Y. 

l=nu 

where  h(i/n)  is  taken  to  be  identically  unity  would  have  a  normal 
limiting  distribution.  However,  (5.4)  is  a  sum  of  independent, 
identically  distributed  random  variables,  and  from  Cram«[l],  Theorem 
23,  we  find  that  a  necessary  condition  that  such  a  sum  have  a  limiting 
normal  distribution  is  the  existance  of  all  moments  of  order  less  than 
two.  Thus,  for  p  >  l/2,  it  is  reasonable  to  say  that  in  general 
S^(p,u,v)  and  S^(p,u,v)  will  not  have  a  limiting  bivariate  normal 
distribution.  What  sort  of  limiting  distribution  these  quantities  do 
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have  when  p  >  l/2,  is  completely  open. 

Thus  from  the  viewpoint  of  being  able  to  say  something  about 
asymptotic  power,  there  is  some  advantage  to  using  the  tests  with 
p  <  1/2  rather  than  those  with  p  >  l/2. 
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6 .  Further  Problems . 

Since  each  memher  of  the  family  S^(p,u,v)  (for  u,v  fixed)  is  a 
consistent  estimator  of  (1/A),  it  follows  that  any  linear  combination 
of  a  finite  number  of  these  estimators  (with  the  weighting  factors 
totaling  unity)  will  also  be  a  consistent  estimator  of  (i/a).  Similarly, 
one  would  gjuess  that  if  H(p)  is  a  probability  distribution  on  [0,l]  then 
.1 

(S' (p,u,v)dH(p) 


I 


should  also  be  a  consistent  estimator  of  (i/A).  It  might  be  worth 
considering  whether  there  is  some  H(p)  which  in  some  sense  gives  a 
"better"  estimator  than  any  individual  S^(p,u,v).  Similar  remarks 
apply  to  S^(p,8).  In  this  case,  the  limiting  behavior  of  S^(p,8) 
as  B  -» 0  would  also  be  of  interest. 

It  is  possible  to  modify  the  two-sample  problem  treated  herein  so 
that  it  becomes  a  two  population  test  of  fit.  Namely,  suppose  that 
H(x)  is  a  given  distribution  function.  Then  under  we  have 

F(x)  =  G(Ax  +  B)  =  H(Cx  +  D)  where  A,  B,  C,  D  are  real  but  un¬ 
specified  constants  (A  >  0,  C  >  O).  The  question  here  is  not  whether 
F(x)  and  G(x)  are  the  same  "type"  of  distribution  but  rather  whether 
they  are  both  the  same  specific  type.  This  modified  problem  can  be 
solved  by  using  the  statistic  proposed  by  Weiss  [J].  let  Z^(x) 

be  defined  for  fixed  u,  v  (0<u<v<l)  by 
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n  X!  [H"^(j/n)3(X  -  X 

(6.1)  Z^(X)  =  - - - 

{  X  h[H"^(j/n)](X.  -  X  )3  ^ 

nu<J<nv  ^  ^ 

and  let  Z^(x)  te  defined  similarly.  Then  from  the  results  of  [5]>  it 
is  easily  seen  that  under  condition  (4.1),  both  of  the  followizig  tests 
are  consistent 

1)  Reject  if  Z^(X)  2^(X)  is  "too  large",  or 

2)  Reject  if  Z^(X)  +  Z^(Y)  is  "too  large". 

Which  of  these  two  tests  is  better  depends  on  their  limiting  power, 
which  in  turn  depends  on  the  unknown  limiting  distributions  of  the 
proposed  statistics. 
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Commander 

San  Antonio  Air  Materiel  Area 
Attn:  Assistant  for  Quality 

Directorate  oF  Maintenance 
Engl:  ering 

Kelly  Air  Fc  :  Baae,  Texas 


CtTmmandcr 

San  Antonio  Air  Materiel  Area 
Attn:  Materiel  QuallU  Division 
Directorate  oF  Supply 
and  Services 
Kelly  AFB,  Texas 


Engineering  Statistics  Group 
Research  Division 
New  York  University 
New  York  53,  hen,  York 


Engineering  Statistics  Unit 
ChemlesF  Corps  Engineering  Agency 


Federal  Telephone  and  Radio  Co. 
100  Kingsland  Road 
Clifton,  New  Jersey 


Commander 

San  Antonio  Air  Materiel  Area 
Attn;  Quality  Control  Division 
Directorate  of  Procurement 
and  Production 

Kelly  APB,  Texas  1 


Field  Inspection  Seetloa 
OES  Division 

Armed  Services  Medical  Procurement 
Agency 

34  Sands  Street 
Brooklyn  1,  New  York 


Logistics  Research  Project 
George  Washington  University 
707  22nd  Street,  N.  W 
Washinglon  7,  D,  C. 


Milwaukee  Air  Procurement  Ditt. 

Altni  Quality  Control  Div. 
770  N.  Plankinton  Avenue 
1  Milwaukee  3,  Wisconsin 


1 


Military  Medical  Stqiply  Agency 
84  Sands  Street 
Brooklyn  1,  New  York 

AUn:  In-Store  Quality  Control 
Materiel  Inspection  Branch 


N.  A.  C.  A. 

IS12  H.  Street,  N.  W. 
Washington  25,  D.  C. 

1  Attnr  Div.  of  Res.  Information 


Commander 

San  Antonio  Air  Materiel  Area 
Attn'  Quality  Control  Office,  SAQ 
Kelly  AF8,  Texas  1 


Commander 

San  Antonio  Air  Materiel  Area 
Attn;  Quality  Control  Division 
Directorate  Special  Weapons 
Kelly  AFB,  Texas  I 


Commander 

Ogden  Air  Materiel  Area 
Aitn;  Assistant  ForQualUy 

Directorate  oF  Maintenance 
Engineering 
Mill  AFG,  Uuh 


1 


Commander 

Ogden  Air  Materiel  Area 

Attn;  Materiel  QuallU  Division 
Directorate  of  Supply 
and  Services 
Hill  AFB,  Uteli 


1 


Mrs.  O.nrothy  Gilford,  Head 

Logistics  and  Mathematical  Statistics  Branch 

Code  436 

Office  of  Naval  Research 

Washlngton25,  0.  C.  3 


Military  Clothing  &  Textile  Supply  Agency 
Philadelphia  Quartermaster  Center 
U.  S .  Amty 

2800  So.  20th  Streel 
Philadelphia  1,  Pa. 

Attn:  Policy  Office  ,  Tedmlcal  Division  1 


Headquarters 

San  Bernardino  Air  MateHel  Area 
U.  S.  Air  Force 
Norton  AFB,  California 

Attn:  Chief,  Planning  6  Control  Offices  1 


Headquarters 

IJ.  S  Army  Signal  Equiprrenl 
Support  Agency 
Fort  Monmouth,  flew  Jersey 

Alin-  $lGFllVES-PP£  1 


Commander 

Ogden  Air  Materiel  Area 

A'tn:  Quality  Control  Division 
Directorate  of  Procurement 
and  Production 
Hill  AFB,  Utah 


Headquarters 

AMC  uillsiic  Missiles  Center,  USAF 
Air  Force  Unit  Post  OfFiee 

Attn*  E,  J.  Lancaster,  Deputy 
for  Quality  Assurance 
Los  Angeles  4S,  California 


Commander 

Ogden  Air  Materiel  Area 

Attn:  Quality  Control  Office,  OOQ 
Hill  AFB,  UUh  2 


David  Taylor  Model  Basin 
Applied  Mathematics  Lab.  (Code  B20) 
V^asIHnglon  7,0.  C. 

Atln;  Dr.  Julius  Lleblein  1 


Detroit  Air  Procurement  Olst . 

Attn:  Quality  Control  Div.  MOKDQ 
W.  Warren  Avenue  &  lonyo  Blvd. 
Oelralt  32/  Michigan 
Attn:  MOHOQI^ 


Director 

National  Security  Agency 
Attn:  REMP-1 
Fort  George  G  Meade 
Maryland 


Director,  Development  Division 
Field  Command 

Armed  Forces  Special  Weapons  Project 
Albuquerque,  New  Mexico 


Eastern  Contract  Management  Region 
Attn:  Directorate,  Quality  Control 
Olmsted  AFB,  Pennsylvania 


Industrial  Division 
Office,  Chief  of  Ordnance 
Department  of  the  Arniy 
Washington  25,  0.  C. 

Attn:  Mr.  Seymour Lorber 


Inspection  and  Quality  Ctxitrol  Div. 
Office,  Asst.  Secretary  ^  Defense 
{$4U 

Washington  25,  0.  C. 

Atln:  Mr.  Irving  B.  Allman 


inspection  and  Quality  C^trol  Div. 
Office,  Aist.  Secret^  Defense 
fSliL) 

Washington  25,  0.  C. 

Atln:  Mr.  JohnJ.  Riordan 


Library 

1  Institute  for  Oefente  Analysis 

Communications  Research  Division 
Von  Neumanrt  Hall 
Princeton,  New  Jersey 


1  Librarian 

Numerical  Analysis  Research 
University  of  California 
Los  Angeles  24,  ^llfornla 

I 

Logistics  Research  Division 
Attn.  MCFR 

Wright-Patterson  AFB,  Ohio 


1 


1 


1 


1 


1 


1 


1 


Newark  Air  Procurement  DIst. 

Attn:  Quality  Control  Division 
218  Market  Street 
Newark ,  Naw  Jersey 


New  York  Air  Prucurement  DIst. 

Attn:  Quality  Control  Div. 
Ill  East  16Ui  SUeet 
New  York  3,  New  York 


Office,  Asst.  Secretary  of  Defense 
fR&E) 

Room  3C1065,  The  Pentagon 
Washington  25,  D.  C. 

Atln:  Technical  Library 


Operations  Analysis  Office 
Headquarters,  Pacific  Air  Forces 
U.  S.  Air  Force,  AP(V<>53 
Attn:  C.  E.  Thompson 
Senior  Analyst 
San  Francisco,  California 


Office,  Chief  of  Engineers 
Department  of  the  Army 
Washinglon  25,  D.  C. 

Attn:  Procurement  Division 
Military  Supply 


Office  of  Opefallcn  Analysis 

DCS/Diseratlons 

Elgin  AFB,  Florida 


Olfice  of  Technical  Services 
Departmtni  of  Commerce 
Washington  25,  0.  C. 


Office  of  the  Chief 
(R&D),  U.  5.  Army 
Arlington  Hall  Station 
Arlington,  Virginia 

Atln:  Dr.  I.  R.  Hershner,  Jr. 


Officer  In  Charge 

U .  S .  Nai^  Central  Torpedo  Office 
Newport,  Rhode  Island 
Attn:  Mr.  C.  B.  Habichl 


Ordnance  Corps 
Industrial  Engineering  Div. 
Diamond  Ordnance  Fuae  Lab. 
Washington  25,  0.  C. 


Ordnance  Mission 
White  Sands  l^ovlng  Ground 
Las  Cruses,  New  Mexico 
Attn:  Mr  Paul  G.  Cox 


Philadelphia  Air  Procurement  Dlit. 

Attn:  QualiW  Control  Division 
1411  Walnut  Street 
Philadelphia  2,  Pennsylvania 
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Physical  Research  Branch 
Evarfs  Signal  Lah.,  SCEL 
Belmar,  New  Jem^y 
Attn:  Mr.  Joseph  Weinstein 


The  RAND  Corporation 
1700  Main  Street 
Santa  Monica,  California 
Attn:  Library 


U.  S.  N.  0.  T.  5. 
Pasioena  Annex 
3202  E.  Foothill  Blvd. 
1  Pasadena  8<  California 
Attn;  Mr.  S.  Caspar 


Western  Contract  Management  Region 
Attn:  Directorate,  Quality  Control 
1  Mira  Loma  AF  Station,  Calilomia 


Dr.  A.  Chames 
The  Technological  Institute 
Northwestern  university 
Evanston,  Illinois 

1 

Mr.  W.  H.  Clatworthy 
Bettis  Plan,  Westinghouse 
Electric  Corporation 
1  Box  1468 

Pittsburgh  30,  Pennsylvania 


Report  Library 
Unlveistty  of  California 
Los  Alamos  Scientific  Lab. 
P.  0.  Box  1663 
Los  Alamos,  New  Mexico 


Rguhel  Oetfelupt'onl  Croup 
R«jStone  Arsenal 
Huntsvill Alabama 

Attn:  ..i.  E.  1..  Bombara 


Or.  Adam  AbruzzI 

Dept,  of  Economics  and  Enalneering 
Stevens  Institute  of  Technology 
Hofxiken,  New  Jersey 


Mr.  William  E.  Gilbert,  Chief 
Maliteriialicat  Statistics  Branch 
Atomic  Energy  Commission 
Washington  25,  D.  C. 


San  Francisco  Air  Procurement  Olst. 

Attn:  QuallW  Control  Division 
Oakland  Army  Termirul,  Bldg.  1 
West  Grand  &  Ma/lllme 
Oakland  14,  California 


Scranton  Or.i()ance  Plant 
156  Ced.’f  Avenue 
Scranton,  Pennsylvania 
Attn:  Mr  Carl  0.  Larsen 
Oil«f  Inspector 


Special  Project  Office 
SP2016 

Department  of  thf  ^'.'vy 
Washington  25,  C .  w. 

At**-  Or.  H.  Weingarien 


$t.indarcls  Branch 
F..t-urement  Division 
DCS  Logistics,  U.  S.  Army 
WarMrglon  25,  0.  C. 

Attn:  Mr.  Silas  Williams, 


Statistical  Lab''  .yy 
University  of  Ca.nornla 
Berkeley  4 .  California 


Technical  Informatlort  Officer 
Naval  Research  Laborat(rr; 
Washington  25.  0 .  C. 


Technical  Operation,  Inc. 

and  C.  0.  R.  G. 

Hg,  Continental  Army  Command 
Fort  Monroe,  Virginia 


U.  S.  Naval  Avionics  Facility 
Indianapolis  18,  irrdiana 
Attn:  Library 


Professor  T.  W.  Arxierson 
Oenartnent  of  Mathematical  Statistics 
Colu«d>la  University 
New  York  27,  New  York 


1 

Professor  Fred  C.  Andrews 
Mathematics  Department 
University  of  Oregon 
Eugene,  Oregon 


1  Dr,  Max  Astraehan 

Department  of  Logistic 
The  RAND  Corporation 
1700  Main  Street 
Santa  Monica,  Callform.i 


Professor  Robert  Bechhefer 
SIMey  School  of  Meet:.  Engineering 
Cornell  University 
Ithaca,  New  York 


,  Or.  R.  E.  Beckwith 

1  Aeronutrortic 

Ford  Road 

Newport  Beach,  California 


1  Professor  J.  N.  BertcllonI 
Western  Reserve  Un'  ^erslly 
Cleveland,  Ohio 


6 


1 


1 


Mr.  Carlton  M.  Beyer 
Office  of  Guided  Missiles 
O'f.ve  of  Asst.  Secretary  of  Defense 
fR&E) 

Washington,  D.  C. 


Professor  Z.  W.  Blmbaum 
'  ^oratory  of  Statistical  Research 
•Htpartmem  of  Mathematics 
University  of  Washington 
Seattle  5,  Wdslilnglon 


U.  S.  Nave'  F.;iglneerlng  Experiment 
Sutton 

Anaoolls,  Maryland 

S^ri:  Mr.  F.  R.  OelPrIare 


U.  S.  Naval  Inspector  of  Ordnance 
400  S.  Beiget  Street 
Mishawaka,  Indiana 


U.  S.  Naval  Inspector  of  Ordnance 
Eastman  Kodak  Company 
50  W.  Main  Street 
Rocliester  14,  New  York 


Professor  Russell  Bradi 
Department  of  Mathematics 
University  of  Kansas 
Lawrence,  Kansas 


Pr^essor  Irving  W.  Burr 
Department  of  Mathematics 
1  Ptmue  University 
Lafayette,  Indiana 


Mr.  G.  Burrows 
Knolls  Atomic  Power  Lab. 
1  Schenectady,  New  York 


Professor  Paul  Clifford 
Nav  Jersey  State  Teachers  College 
1  Montclaire,  New  Jersey 


Professor  W.  G .  Cochran 
Department  of  Statistics 
Harvard  University 
1  Cambridge,  Massachusetts 


Professor  C.  C.  Cockerham 
Institute  of  Statistics 
State  College  Section 
1  North  Carolina  State  College 
Raleigh,  Ncrih  Carolina 


Professor  Edward  P.  Coleman 
Engineering  Department 
1  University  of  California 

Los  Angeles  24,  California 


Dr.  Louis  Court 
Division  17 

National  Bureau  of  Standards 
1  Washington  25,  D.  C. 


Miss  Gertrude  M.  Cox 
Instllute  of  Statistles 
North  Carol Irta  Slate  College 
1  Raleigh,  North  Carolina 


1 


1 


Or.  Joseph  Daly 
U.  S.  Census  Bureau 
Washington  25,  D.  C. 


Professor  Cyrus  Derman 
Department  of  Industrial  EngInceiTng 
Columbia  University 
New  York  27,  New  York 


Mr.  M.  F.  Dodge 
Rutgers  University 
New  Brunswick ,  New  Jersey 


1 

Dr.  FrariCls  Dreseh 
Stanford  Research  Institute 
3.13  Ravenswood  Avenue 
Menlo  Park,  California 


1  Professor  Aeheson  J.  Duncan 

Department  of  Industrial  Engineering 
Johns  Hopkins  University 
Baltimore  18,  Maryland 


1  Professor  Mever  Dwass 
Depertment  of  Mathemittci 
Northwestern  University 
Evanston,  Illinois 


1 


I 


ProleMorD.A.S.  Fraw 
Deparunent  of  M*Ui.ntttlct 
Unlv.r»IW  o!  Twonto 
Toronto  5,  Canada 


1 


1 


1 


1 
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1 


1 


1 


1 


1 


1 


1 


1 


1 


U.  S.  Naval  Inspector  of  Ordnance 
Statistics  Department 
Naval  Gun  Faelora 
Washington,  0.  C. 

Attn:  C.  0.  Hock 


C^alnA.  E.  Chapman 

BMing  Aircraft  Corporation 
Seattle,  Washington 
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Mr.  Beinard  P.  Goldsmith 
AsiMiate  Professor 
Northwestern  University 
Huntl[t9lon  Avenue 
Boston  IS ,  Massachusetts 


Professor  Leo  A.  Goodman 
Statistical  Research  Center 
University  of  Chicago 
Chicago  i7,  Illinois 


Professor  W.  G.  Jones 
Department  of  Stattstles 
Michigan  State  University 
East  Lansing*  Michigan 


Mr.  J.  P.  Kearney 
Quality  Control  Division 
General  Services  Administration 
Room  6316*  Region  3  Building 
1  Washington  25*  D.  C. 


D.  E .  Newnham 

Chief*  Industrial  Engineering  Division 
Comptroller 

1  Hd<M .  *  San  Bernardino  Air  Materiel 
Area 

Norton  Air  Force  Base*  California  1 


1 


Professor  J.  Ncyman 
Department  of  Statistics 
University  of  California 
Berkeley  4,  Califwnia 


1 


Mr.  Leon  Gilford 
Operations  Research  Inc. 
3605  Cameron  Street 
Silver  Springs*  Maryland 


Professor  Oscar  Kempthome 
Statistics  Laboratory 
Iowa  State  College 
Ames  *  Iowa 


Or .  J .  Greenwood 
Directorate  of  Intelligence 
Hq.*  U  S.  Air  Force 
Washington  25*  D.  C. 


Professor  Solomon  Kullhack 
Deportment  of  Statistics 
Gcorje  Washington  University 
1  WashirAon  7*  D.  C. 


Professor  Frank  M.  Gryna*  Jr. 
University  College 
Rutgers  University 
New  Brunswick,  New  Jersey 


1 


Dr.  Carl  F.  KoiSKk 
Statistlos  and  Operations  Research 
Research  Cerrtcr  -  IBM 
P.  0.  Box  Hfi*  Lamb  Estate 
Yorktown  Heights*  New  York 


Dr.  Donald  Guthrie 
Stanford  Research  Institute 
333  Ravenswood  Avenue 
Menlo  Park*  California 


1 


Mr.  Howard  Laltin 
The  RAND  Corporation 
1700  Main  Street 
Sarda  Monica,  Caflfomta 


Or.  Theodore  E.  Harris 
The  RAND  Cotvoralion 
1700  Main  Street 
Santa  Monica,  California 


Or.  E.  L.  LeCferv*  Chief 
Biometrical  Services 

1  Agrlcultieal  Research  Service*  USDA 
BetUvIlie*  Maryland 


Or.  Leon  H.  Herhaek 
Department  of  Mathematics 
New  York  Univenily 
New  York  3,  New  York 


Professor  Sebastian  B.  LlUeuer 
411  Engineering  Building 
1  Columbia  University 

New  York  27*  New  York 


Professor  W.  HIrsch 

Institute  of  Mathematloal  Sciences 

New  York  University 

New  York  3*  New  York 


Or.  Eugene  Lukact 
Department  of  Mathematics 
1  Catholic  University 

Washington  17*  0.  C. 


Or .  Paul  G .  HocI 
Department  of  MathematUi 
University  of  California 
Los  Angeles  24,  California 


Professor  Harold  Hotelling 
Associate  Director 
Institute  of  Statistics 
University  of  North  Carolina 
Chapel  Hill,  North  Carolina 


Dr.  Rot>ert  Lundtgard 
Logistics  and  Mathematical 
Statistics  franch 
Office  of  Naval  Research 
Washington  25,  D.  C- 


Professor  Frank  Massey 
Schod  of  Public  Health 
University  of  California 
Los  Angeles  24,  California 


Professor  L.  Hiewlcr 
School  of  Business  Administration 
University  of  Minnesota 
Minneapolis*  Minnesota 


Professor  G.  W.  McElrath 
Department  of  industrial  Engineering 
Univinily  of  Minnesota 
Minneapolis  14*  Minnesota 


Mr.  Rudolf  Husser 
Visiting  Research  MatNmstlclant 
Numerical  Analysis  Research 
Univenily  of  ullfomia 
Los  Angeles  24*  California 


Dr.  Paul  Meyer 
Department  of  Mathematics 
Washington  State  University 
Pullman,  Washli^on 


Dr.  James  R.  Jackson 
Management  Sciences  Research  Project 
63  Aemli  iitratlorr  &lldlr^ 

University  of  California 
Los  Angeles  24*  California 


Dr.  W.  C.  Jacob 
Agronomy  Department 
Universify  of  Illinois 
Urbana*  Illinois 


Dear)  Paul  E.  Mohn 
School  of  EnglnecrNr^ 
Ufllvffsl^  of  Buffalo 
Buffalo*  New  York 


Mr.  R.  S.  Murphy 

Bel!  Telephone  Laboratories*  Inc. 

463  West  Street 

New  York  14*  New  York 


1 


Mr.  Monroe  Norden 
Research  Division 
College  of  Engineering 
New  York  University 
New  York  53*  New  York 


1 


1 


Mr .  Fred  Okano 

National  Aeronautics  A  Space  AdmInIsUatlon 
Reliability  It  Systems  Analysis  Office 
Room  A-131 
1520  “H"  Street,  N.  W. 

Washington  25*  D.  C. 


1 


1  Professor  £.  G.  Olds 

Department  of  Mathematics 

Carnegie  Inslllute  of  Technology 

Plltsb^h  13*  Pennsylvania  1 


1  Professor  £.  R.  Otl 

Department  of  Mathematics 
Rutgers  University 

New  Brunswick,  New  Jersey  1 


1  Mr.  Cyril  Peckham 
Project  Globe 
University  of  Dayton 

Dayton*  Ohio  1 


1  Dr.  Richard  Post 

Department  of  Mathematics 
San  Jest  State  College 
San  Jose,  California 


1  ProfesserP.  H.  Randolph 
Purdue  University 

Department  of  Industrial  Engineering 
Lafayelle,  Indiana  1 


1  Professor  George  J  Resnikoff 

Oepartiiwnt  of  Industrial  Engineering 

Illinois  Institute  of  Technology 

Chicago  16*  Illinois  1 


1  Dr.  Paul  R.  Rider 
Chief  Statistician 

Aeronautical  Research  Lab.*  WADC 
Wright-Patterson  AFB,  Ohio  1 


1  Professor  Herbert  Robbins 

Mathematical  Statistics  Department 
Columbia  University 

New  York  27 ,  New  York  1 


1 


Df.  Harry  G.  Romlg 
351  Alma  Real  Drive 
Pacific  Palisades*  California 


1 


I 


Dr.  Harry  Rosenblatt 
Statistical  Research  OIvIsIoa 
Bureau  of  Census 
Washington  25*  0.  C. 


1 


Professor  Murray  Rosenblntl 
1  Department  of  Mathematics 
Brown  University 

Providence,  Rhode  Island  1 
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Professor  Herman  Robin 
Department  of  Statistics 
Michigan  State  University 
East  Lansln9,  Michigan 


Professor  Norman  Rudy 
Statistics  Department 
Sacramento  State  College 
Sacraments,  California 


Miss  Marlon  M.  S^ndomlre 
U.  3.  Department  of  Agriculture 
Western  Regional  Laboratory 
Biometrical  Services 
Albany  10,  California 


Professor  I.  Richard  Savage 
School  of  Business  Administration 
University  of  Minnesota 
Minneapolis,  Minnesota 


Professor  L.  J.  Savage 
Mathematirs  Department 
University  of  Michigan 
Ann  Arbor,  Michigan 


Professor  Henrv  Scheff/ 
Department  of  Statistics 
University  of  California 
BerVe-fey  A,  Calirnmla 


Professor  Robert  Sehlalfer 
Graduate  School  of  Business 
Administration 
Harvard  University 
Boston  63,  Massachusetts 


Professor  Seymour  Sherman 
Dep-vtment  of  Maihematlcs 
Wayne  Slate  University 
Oeuolt  2,  Michigan 


Mr.  Walter  Shewhart 

Bell  Telephone  Laboratories,  Inc. 

Murray  Hill,  New  Jersey 


Or.  Rosedith  Sltgreaves 
Teacheri  College 
Columbia  University 
New  York  27,  Mew  York 


Dr.  Milton  Sobel 
Statistics  Department 
University  of  Minnesota 
Minneapolis,  Minnesota 


Professor  Prank  Spilaer 
Department  of  Mathematics 
University  of  Minnesota 
Minrseapolls,  Minnesota 


Mr.  Arthur  Stetn 
Cornell  Aerorautkai  Lab.,  Inc. 
P.  0.  Bom  239 
Buffalo  21,  Now  York 


Professor  W.  Allen  Wallis 
School  of  Business 
University  ol  Chicago 
Chicago 37,  Illinois 


Dr.  Irving  Weiss 
The  Mitre  Corporation 
Bedford,  Massaehuitlts 


Captain  Burton  L.  Welter 
AFPR  Office 

Miirtln  Aircraft  Corporation 
Denver,  Colorado 


Professor  Mason  Wescott 

Editor,  Industrial  Quality  Control 

Other  Porelon  Addresses 

Rutgers  Univertitv 

1 

New  Brunswick,  New  Jersey 

1 

Professor  Maurice  H.  Bella 

University  of  Melbourne 

Professors.  S.  Wilks 

Department  of  Mathematics 

Princeton  University 

Carlton  N.  3 

Victoria,  Australia 

1 

Princeton,  New  Jersey 

1 

Professor  Toslo  Kitagawa 

Math^tlcal  Institute 

Professor  J.  Wolfowllc 

Department  of  Mathematics 

Cornell  University 

Faculty  of  Sclarct 

Kyusyu  University 

Fukuoka,  Japan 

Ithaca,  New  York 

1 

1 

Kanlehl  KMonagl 

1 

Dr.  Max  A.  Woodbury 

Oeptftment  of  Mathematics 

College  ^  Engineering 

New  York  University 

New  Ved  S3,  New  York 

1 

Managing  Director 

UnlArt  of  Japarme  Scientists 

4t  Engineers 

2,  1-chome,  Cyobaihl 

Chu9*iii>,  Tokyo,  Japan 

Dr.  Lai  Vtrman 

1 

Distribution  via  ONR  London 

Director 

Indian  Standards  Inslltute 

New  Delhi,  India 

Commanding  Officer 

Branch  Offtet 

Navy  No.  100 

Fleet  Post  Office 

New  York.  New  York 

ProfeikorP.  C.  Mahatanobis 

1 

Indian  Stallitleal  Instllule 

20)  Oarrackpore  Trunk  Read 

Calcutta  35,  India 

Or.  WiiMan  R.  Buckland 

22  Ryder  Sttret 

London  S.  W,  1 
£r.^and 

N.  T.  Mathew,  SRC  Unit 

India  Slatlitieal  Institute 

1 

8  King  George  Road 

Niw  Delhi,  fi^ii 

1 

Professor  Georges  Darmols 

OlfKtor,  Instllule  de  Stailstlgue 
Unlvrrsity  «f  Paris 

11  Rue  rltfff  Curie 

ProfHiorSigelil  MerigMI 

Syean  Mlnamlmall  6 

Suolnunl'liu,  Tokyo,  Japan 

1 

Paris  S,  France 

1 

R.  G.  Narailmhan 

1 

Professor  R.  Fertet 

Instltut  Harei  Poincare 

Paris ,  Frar^e 

1 

SRC  Unit  Calevtia 

Indian  Siaiiiiloai  Instttuie 

9B  Csptanado  Cast 

Calcutta  1,  liwla 

1 

Dr.  Geoffrey  Gregory 

4,  Osborne  Grove 
tiatliv,  Lheadic 

Chesnire,  England 

1 

H.  D.  Shourle 

Dirorinr 

National  Preduellvlly  Ceunoll 
(lelf  Links 

A.  Held 

New  Dolhl,  India 

1 

rcrcsovel  83 

Viruffl,  Denmark 

1 

P,  V.  Slvaramkrlihan 

SQC  Unll,  Indian  flatlilloal  Inst 

QuHns  Road,  Oevernnent  Oflleei 

Dr.  H.  C.  Hamakar 

Philips  Research  LehoratcrFei 
Eideidiovtn,  Ntiherland', 

1 

ftrliding 

Bombay,  India 

Srinagabhusharti 

SQUUnit  toigalivo 
^SKfJ  THHnologleal  hti. 

Bangalorf  1,  India 

1 

Mr.  1.  0.  Hlli 

SlalUtleai  Advisory  Unit 

Mlnisinr  of  Supply 

Lo^on  W.  1.  England 

1 

Pruftsso*  CdusHe  Valanivela 

1 

f^essor  M.  0.  Kendall 

London  School  of  Economies 

London,  England 

1 

Rtpubliea  917-CIEF 

Santiago,  Chile 

Mr.  Cosmo  Villegas 

X 

Profosior  A.  Wallhor 

Teehniicho  Hochishulo 

Darmstadt,  Germany 

1 

Inst,  tie  Matemallea  y  Esiadlstloa 

Av.  J,  Herrera  y  ReitiHI 

Montevideo,  Uruguay 

Additional  cepits  far  preiKt  laadars 
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